Optimal universal quantum cloning and state estimation 
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We derive a tight upper bound for the fidehty of a universal N —> AI qubit doner, valid for any 

M > A'^, where the output of the doner is required to be supported on the symmetric subspace. Our 

proof is based on the concatenation of two doners and the connection between quantum cloning 

. . . and quantum state estimation. We generalise the operation of a quantum doner to mixed and/or 

00 ' entangled input qubits described by a density matrix supported on the symmetric subspace of the 

O^ I constituent qubits. We also extend the validity of optimal state estimation methods to inputs of 

0\ , this kind. 

^-H , 
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^^ Perfect quantum cloning is impossible W. This notwithstanding, we may ask how well we can clone quantum states. 

_— . Buzek and Hillery, who were the first to address this problem, provided an example of a quantum device which can 

clone an unknown pure state of a single input qubit (a two-state system) into two output qubits, albeit with a certain 
^f^ fidelity smaller than one |Q. Their construction was subsequently shown to be optimal by Brufi et al. [^. In this 
^ letter we derive the optimal fidelity of a universal and symmetric quantum cloning machine (QCM) which acts on N 
0^ ' original qubits and generates M clones. 

A universal N —^ M quantum doner is a quantum machine which performs a prescribed unitary transformation on 

an extended input which contains N original qubits, M — N "blank" qubits and K auxiliary qubits, and which outputs 

M clones together with the K auxiliary qubits. The original qubits are all in the same (unknown and pure) quantum 

h — , state described by the density operator g"^ = i(]l + s*" • (?), where s™ is the original Bloch vector. Both "blanks" 

^\ •■ and the auxiliary qubits are initially in some prescribed quantum state. The output qubits are in an entangled state 

and in the present work we require that the density operator describing the state of the M clones is supported on the 

symmetric subspace. This guarantees that all the output qubits are indistinguishable and in the same state described 

by the reduced density operator p°"*. We comment on relaxing this assumption at the end of the paper. 

S It has been shown that universal 1 — > 2 doners can only shrink the original Bloch vector, without changing 

C^ its orientation in the Bloch sphere 0. The same argument as given in 0] (namely the impossibility to find a 

^ , transformation that rotates any Bloch vector of the 1-particle reduced density matrix by the same angle) applies also 

<J ■■ generally for an iV ^ M qubit doner. Therefore, the operation of a universal QCM can be characterised by the 

;> ] shrinking factor ri{N, M), which is also known in the literature as the Black Cow factor [H, and the reduced output 

density operator is of the form p°"* = ^(1 + rj{N, M)s™ ■ a). Universal N ^ M quantum cloning machines may be 

constructed in many diff'erent ways, the best constructions are those which maximize ri{N, M) (i.e. which maximize 

^ . the fidelity of the cloning machine) and we refer to them as the optimal doners. 

■ - ' Gisin and Massar have constructed a class of universal N ^ M QCMs and showed, using numerical methods, that 

for A^ < 7 their doners are optimal |g]. Our derivation of the upper bound for ri{N, M) is quite general and does not 

refer to any specific realisation of the universal cloning machines. In particular it shows that the Gisin-Massar doners 

saturate this bound for any TV and M > N . Our approach avoids an elaborate optimisation procedure, extends the 

class of allowed inputs to mixed and/or entangled states of the original qubits which belong to the symmetric subspace 

and sheds some light on the connection between optimal quantum cloning and optimal quantum state estimation. 

The proof is based on the concatenation of two quantum doners and on associating the upper bound on the fidelity 

of an 7\f ^ cx) doner with the fidelity of the optimal state estimation of M qubits, given in |^. 

We concatenate two cloning machines in the following way. The first doner is an TV ^ M universal machine 
characterised by the shrinking factor rj{N,M). The M clones from the output of the first doner are then taken as 
originals for the input into the second cloning machine which creates infinitely many clones with the shrinking factor 
ri{M, cxd). We now write down two statements which will be proved after unfolding the main result: 

a) The shrinking factors for concatenated doners multiply. 

b) The equality 

r^"^l,^,{M,^)=nTeas{M) (1) 
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holds. Here tj'^I^^^M) corresponds to the optimal state estimation derived in Q, and its meaning will be 
explained below. 

Due to statement a), the shrinking factors of universal cloning machines in sequence multiply. Moreover, the 
sequence of the two machines cannot perform better than the optimal iV — > oo universal doner, otherwise the iV — > cx) 
universal doner would not be optimal. Thus we arrive at the following inequality: 

i1qcm{N,M)-t^qcm{M,^)<ti"^'cm{N,^) . (2) 

This means that the lowest upper bound for the general N ^f M doner is given by 



Vqcm{N,M)< '1--Z: [ (3) 



We have thus reduced the optimality problem of the N ^ M doner to the task of finding the optimal N ^ oo doner. 
Now we can use statement b) and the explicit form of fj'^easi^^) (see 0), namely 

vTeasiM) = J^ (4) 



to conclude the central result that for any M > N 



N M- 



^ZcMiN,M)^-'-^. (5) 

For pure input states this corresponds to the optimal fidelity 

which is achieved by the cloning transformations proposed in ||^ . (For g*" ~ \ ■0} ("0 I the fidelity is defined as 

Let us note in passing that as the consequence of the factorisation property (J^) we can produce M clones from N 
originals either by applying directly the optimal N —> M doner or by taking any number of intermediate steps in 
order to realise the cloning process, using the optimal transformation at each step; both ways lead to the same overall 
shrinking factor. 

Let us now justify statements a) and b). 

In order to prove a) we describe an N -^ M doner in terms of a completely positive map C nm which maps input 
density operators of N identical pure originals into output density operators of M clones, such that for any state 
|'0)(''/'l of 8- single input qubit (original) we have 

TrM-i[CNMm{M'''')]^^{N,Mm{^\ + {l-r^{N,M))h, (7) 

where the trace is performed on any M — 1 qubits (for an overview of completely positive operators see W). 

Let qn be a density operator of N qubits which is supported on the symmetric subspace of the 2^ dimensional 
Hilbert space. We can always write gjv as a linear combination of direct products of identical pure states, gjv — 
Yl,i^i\'^i){'^i\'^^ , where Yli'^i = 1; ^-^^ we do not require that all values a^ are positive [§J|]. The linearity of the 
completely positive map and its universality, i.e. the fact that rj{N,M) does not depend on \il}), allow us to extend 
Eq.(^ to the more general form 

TrM-i[CNM{QN)] ^ v{N,M)q+ {I - ri{N,M))U. (8) 

where g = Trjv-i[£'Ar]. Now, suppose we concatenate an iV ^ M and an M — » L doner and view it as an TV — > L 
doner. It evolves the initial N qubit state qn first into the M qubit state g^ and then into the L qubit state p^. The 
corresponding single qubit reduced density operators are g, g ^ T^'<^m-i[Qm] ^^"^ 6' = Tri_i[p^] . Following Eq.(||) 
we can write 

/ = v{M,L)g + (1 - viM,L))^l =ri{N,M) ■ rj{M,L)g+{l - v{N,M) ■ v{M,L))^l (9) 



i.e. indeed r]{N, L) = r]{N, M) ■ r]{M, L). 

We will now prove statement b). Equation (H) was obtained assuming the following result (due to B): given M 
qubits all in an unknown quantum state | f/;) there exists a universal POVM measurement {P^} [^ which leads to 
the best possible estimation of | ip) with fidelity F{M) = Ifx^, or, equivalently, with fi{M) — -^r^ ■ The outcome of 

each instance of the measurement provides, with probability Pf_i{ip) — Tr(P^|i/))(V'|®^), the "candidate" | -f/'/i) for | ^). 
The fidelity Fmeas{M) is then calculated from the outcomes of the measurement as 

FmeasiM) = Y,PM)Mi'^)\^ = (^1^1 ^), (10) 

where g = X^„P/j(V')I''/'m)(''/'mI ■ l'^ ^^^ optimal, universal state estimating procedure the fidelity must not depend on 
"0, thus g can also be written as 

= fiTeas{Mm{ij\ + (1 - C^L(M))^1- (11) 

The optimal measurement of this type can be viewed as an M —>■ cxj doner because after reading each outcome 
we can prepare any number of "candidates" , in particular infinitely many of them, with the average reconstruction 
fidelity F^^^g{M) with respect to the originals. Clearly this procedure cannot provide a larger shrinking factor than 
the optimal M —>■ L doner and we find 

riTeas{M)<vTcM{M.L) (12) 

for any L > Af , in particular for L —>■ oo. 

Let us now show that for L ^ oo the formula (O) becomes the equality. To see this let us concatenate an M ^ i 
doner with a subsequent optimal state estimating measurement. The input to the doner is of the form |?/j)('(/'|®^^ 
and the output is described by the density operator gL which is of the form J2i cii\ipi){''Pi\'^^ , where J2i cti = 1. The 
reduced density operator of each output qubit is g = Tr^-igL = J2i^i\''Pi){^i\ — '?(-^: ^)l^)(^l + (1 ~ vi^iL,))^!. 
The doner M ^^ L concatenated with the state estimation on L qubits can be viewed as the state estimation performed 
on M qubits. The total procedure gives the fidelity of estimating Jt/i) which can be written as 

FmeasiM) = (^| ^ Tr(P^ei) |^^) (0^ J V') = J] (V'|a.Tr(Pj^,) (^, |®^) |^^) (V' J V') (13) 

= E(^l"^[^-as(^)IV^.)('A.| + (1 - fjTeas{L))lim (14) 

i 

which for i ^ oo becomes (due to VmlasiF) -^ 1) 

Pmea.(M) ^E(V'|a.|V'^)(0^|V') = (V^lffl^) = -{I + f]QC m{M , (x)). (15) 

i 

The concatenation of a doner with a measurement cannot perform better than the optimal measurement, thus we 
can write 

l^f^^{M,^)<fiTeas{M) . (16) 

Combining equations ( |l^ ) and (16) finally leads to 

77^P^,,(Af,^)=Cfe*a.(M) , (17) 

thus proving statement b). 

Before concluding, we want to stress that, as a consequence of what was shown above, we can extend the operation 
of a cloning machine and the application of an optimal measurement to any input density operator of N qubits which 
has support on the symmetric subspace. The properties of the universal cloning machine as defined at the beginning 
of this paper allow us to conclude that the same machine can operate on any such symmetric density operator and 
shrinks the Bloch vector of the reduced input density matrix by a fixed amount, independent of the initial length. 
Notice also that the optimal machine, for products of pure inputs specified by the shrinking factor (ra), is still optimal 
for this extended class of inputs. Actually, if a better cloning machine existed for mixed input states, we would use it 



as the second doner M ^ oo in Eq. (|^), giving a smaller lower bound in Eq. (^). This would lead to a contradiction 
because we already know that universal doners for pure states saturating the bound (|5[) exist g . 

One may want to relax our restriction and consider quantum doners which produce identical clones (i.e. with the 
same single-qubit density operator p°"*), but for which the state of all outputs does not belong to the symmetric 
subspace. This case, in principle, may provide a higher shrinking factor, however, we could neither prove nor disprove 
this with our approach. We leave this problem as a challenge to other colleagues. 

In a similar way as for the doner we can extend the validity of an optimal universal measurement procedure to 
inputs from the symmetric subspace. In this case the goal is to find the optimal measurement for the reduced density 
operator for each input copy. Since we require the process to be universal, we know that the reduced density operator 



reconstructed as the result of the measurement given in Eq. (11) is just the shrunk version of the initial one. We can 
then describe the quality of the procedure in terms of the shrinking factor. We conclude that the optimal measurement 
derived in Rcf. Iq] is also optimal for any input symmetric state. Actually, if this were not the case, we could devise a 
measurement procedure on N initial pure qubits by first applying an iV — > M doner and then an optimal measurement 
on the mixed state of the output M clones. If this global measurement were better than the optimal one of Ref. g] 
we would then obtain with the above procedure a universal measurement for pure states which performs better than 
the one in Ref. [pj, thus finding a contradiction. 

Let us mention in passing that in our discussion we found it convenient to use the shrinking factor, because it has 
an intuitive geometrical meaning both for pure and mixed states, however, one can rephrase the optimality argument 
for universal operations using, for example, the Uhlmann fidelity |10[| for the reduced density operators. 

In conclusion, we have derived the optimal shrinking factor/fidelity for a universal N —^ M doner and generalised 
its operation to a more general case of mixed and/or entangled input states which belong to the symmetric subspace. 
Furthermore we have established the connection between optimal quantum state estimation and optimal quantum 
cloning which allowed us to extend the validity of the optimal state estimation methods B to inputs of the above 
form. 
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